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Abstract 

Quantum conductance calculations on the mechanically deformed monolayers of 
M0S2 and WS2 were performed using the non-equlibrium Green's functions method 
combined with the Landauer-Biittiker approach for ballistic transport together with 
the density-functional based tight binding (DFTB) method. Tensile strain and com- 
pression causes significant changes in the electronic structure of TMD single layers and 
eventually the transition semiconductor-metal occurs for elongations as large as 11% 
for the 2D-isotropic deformations in the hexagonal structure. This transition enhances 
the electron transport in otherwise semiconducting materials. 

1 Introduction 

The semiconductor industry rapidly approaches the performance limit of current silicon- 
based technology. Short-channel effects and gate leakage increase levels of heat dissipation 
as transistors are miniaturized. Silicon, being a 3D material, suffers stability issues and 
quantum size effects if layers approach a critical thinness of a few A. Possible long-term 
solutions to this problem involve replacing silicon with another material that would perform 
better at smaller scales. 2D materials are very attractive alternatives, and the recent progress 
in the field of graphene has demonstrated that nanotechnology is currently in the position 
to process 2D materials. The move to 2D materials also enables new functionalities, such as 
surface-coated devices or flexible electronics. 

The drawback of graphene as a silicon alternative is the absence of a band gap and 
the difficulty to introduce one. There is, however, a rich variety of semiconducting 2D 




Figure 1: (Online color) The atomic structure (a) of 2D transition- metal dichalcogenides of 
TX 2 type in the rectangular representation. The dashed rectangles indicate the units in a 
and b directions. The hexagonal unit cell is highlighted. Schematic representation of the 
electronic transport in the (b) zigzag and (c) armchair directions shown together with the 
possible mechanical deformations: 2D-isotropic (e), uni-axial and parallel to the transport 
axis (ey), uni-axial and perpendicular to the transport axis (ej_). The electrodes, L and R, 
are highlighted with shaded rectangles). 

materials, most prominent are transition-metal dichalcogenides (TMDs) of TX 2 type (where 
T - transition-metal atom, X - chalcogen atom). They are composed of two-dimensional 
X-T-X sheets stacked on top of one another, as shown in Fig. 1. Each T atom is covalently 
bonded to six X atoms located in the top and bottom layers forming a sandwiched material. 
In the bulk, adjacent slabs are held together by weak interlayer bonds, what leads to a quasi- 
2D character of TX 2 . Exfoliation of the bulk down to the monolayer, however, changes the 
electronic properties of the materials by quantum confinement. 1-3 TMD monolayers possess 
high mobility of charge carriers, high thermal stability, and the presence of a distinct band 
gap. These intrinsic properties make them attractive for electronic applications. 

Recently, Kis and coworkers demonstrated that these materials are indeed suitable for 
producing nanoelectrodevices, and the first field-effect transistors, logical circuits and am- 
plifiers have been manufactured. 4-7 Large-scale production appears to be feasible, as liquid 
exfoliation of TMDs is able to produce large area two-dimensional monolayers. 8 

Bulk layered TMDs can be either semiconductors (e.g. T=Mo, W) or metals (e.g. T=Nb, 
Re) , depending on their stoichiometry. It has recently been shown that the former exfoliated 
TMDs in the form of monolayers undergo a transition of their electronic structure. 1,2 While 
bulk structures are indirect band gap semiconductors, the 2D monolayers are direct band 
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gap materials with significantly larger energy gaps of about 2 eV. 2,9 The electronic struc- 
ture of these TMDs can be further tuned by means of nanotube formation, 10 substitutional 
doping 11 ' 12 or mechanical deformations. 13-17 The latter one is particularly interesting as the 
semiconductor- metal transition occurs at elongations of about 11% without bond breaking. 
Such rather large elongations have been observed in the experiment, e.g. in tensile tests of 
WS 2 nanotubes. 18 

The proposed semiconductor- metal transition 13-17 has been discussed so far solely by 
band structure calculations. However, stretched semiconductors are likely to expose states 
at the Fermi level due to broken (dangling) bonds, or bonds that are about to break, that 
do not confirm electrical conductance. It is imperative to study the electron conductance 
through these materials in order to conclude if 2D TS2 layers could find applications as 
nanoelectromechanical switches. Here, we analyze strain effects on the electronic transport 
properties of two-dimensional TMDs using density functional based methods and 2D quan- 
tum transport calculations. 

2 Methods 

The fully optimized structures of TS 2 with T = Mo, W, and Nb were taken from our 
previous work. 2 For computational details and optimized lattice parameters of TMDs, the 
reader is referred to the work of Kuc et al. 2 In short: all stability and band structure 
calculations have been carried out using density-functional theory (DFT) employing the 
PBE (Perdew-Burke-Ernzerhof) functional. 19 The equilibrium structures were exposed to 
mechanical strain following three possible strategies: 1) 2D-isotropic deformation (e) on 
the hexagonal unit cells with simultaneous change of a and b lattice vectors, 2) uni-axial 
deformation of a along the armchair edges (e n ) in rectangular unit cell, and 3) uni-axial 
deformation of b along the zigzag edges (e&) in the rectangular unit cell. Depending on 
the direction of the current, the uni-axial deformations can be either perpendicular (e±) or 
parallel (e||). For example, the notation of e a ± denotes the uni-axial strain along a lattice 
vector and perpendicular to Q (see Fig. 1). Note, that the uni-axial tensile strain along the 
b and a vectors correspond also to the deformations of large-diameter zigzag and armchair 
TMD nanotubes, respectively. The strain is defined as e = L ~^ , where L and L are 
equilibrium and strained lattice values. 

The coherent electronic transport calculations were carried out using our in-house non- 
equilibrium Green's functions (NEGF) 20 code combined with the Landauer-Buttiker ap- 
proach, that takes advantage of the density functional based tight-binding (DFTB) Hamil- 
tonians. 21-23 The Slater-Koster parameters for DFTB calculations were created following 
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the approach described by Oliveiara et al. 22 and Seifert et al., 23 and have been validated by 
band structure comparison with DFT calculations. [For the band structures calculated at the 
DFT/PBE level see Supporting Information]. 24 We have used the PBE 19 density- functional, 
numerical basis functions, and the quadratic confinement potential with r = 5.5, 5.8, 5.9, 
and 3.8 Bohr for Mo, W, Nb, and S, respectively. DFTB band structure were calculated 
using DFTB+ code. 25 The mean free path of the carriers in TMDs is greater than the sam- 
ple dimensions and their high mobility 4 implies that the carriers can move long distances 
without scattering. Hence, it is expected that the ballistic approximation is appropriate for 
studying transport in TMDs. 26 

The electron conductivity was calculated using a system that consists of a finite TMD 
monolayer as scattering region (conductor), which is connected to two semi- infinite TMD 
monolayers that act as electrodes, R and L (see Fig. I). In this way, the so-called bulk 
conductivity is simulated, where the interface between the electrodes and the conductor is 
well defined. In the bulk conductivity, the scattering region has to be large enough to avoid 
the interaction between R and L. We have converged the scattering region to seven unit 
cells between the electrodes. The direction perpendicular to the transport axis is assumed 
to be two dimensional by accounting that each border atom in the unit cell is surrounded by 
translational equivalence of atoms according to the symmetry of the periodic system. In this 
direction, the periodicity is described by the T-point approximation, recently justified for 
analogues systems. 27 To avoid spurious overlap between atoms from the image cells, we have 
used twelve unit cells in the direction perpendicular to the current. In this way, the edge 
effects arising from the typical one-dimensional representation of the layers as nanoribbons 
are discarded from the simulations. At low temperature, the electron conductivity through 
such a system can be described using Fisher-Lee Formula: 28 



the coupling matrices. Self-energies, Sl,_r, which describe the influence of the semi-infinite 
electrodes on the electronic structure of the device, can be calculated via a self-consistent 
procedure. 29 

Hirshfeld charge calculations were performed using the ADF-BAND code, with PBE 
functional and triple-zeta basis sets with one polarization function (TZP), including scalar 
ZORA relativistic correction. 30 

The elastic properties of the monolayers under the uni-axial deformations correspond to 
those of TMD nanotubes with large diameters. Therefore, we have used the tensile stress as 
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force F acting on the area A. For the area calculation, we have used the fixed lattice vector 
(e.g. if a was deformed, b was used to calculate A), multiplied with the optimized monolayer 
thickness of 6.2A. 

3 Results and Discussion 

Tensile stress and compression within the monolayers cause changes in the geometry, 
namely the elongation or shrinking of the bond lengths and bond angles. TMDs can be 
exposed to such strong mechanical deformations without breaking the T-S bonds due to 
very flexible S-T-S bond angles, which can be easily stretched or compressed by up to 
about 10°. Figs. 2 a and b show the calculated interatomic distances (T-S) and bond angles 
(S-T-S) as function of the applied strain. T-S bond distances become longer by about 
0.01 A per 1.5% isotropic tensile strain. At the same time, significant changes are obtained 
in the S-T-S bond angles: about 1.5° per 1.5% isotropic tensile strain. Both, the bond 
lengths and angles scale nearly linearly with the applied deformation for up to about 11% 
strain/compression. Uni-axial deformation cause slower changes in the geometry for a given 
e comparing with the isotropic deformations. 

Changes in the geometry of TMD monolayers under mechanical deformations also affect 
the electronic structure of these materials. For isotropic deformations, we have calculated 
evolution of the atomic charge (q) for the T atoms using Hirshfeld atomic charge analysis. 
The Mo and W atomic charges are shown in Fig. 2 c as function of e. The transition 
metal atoms become more positive with increasing the tensile strain, at the same time the 
sulfur atoms get more negative. This leads to a charge depletion around the T atoms and 
consequently charge excess around S atoms. Atomic charge analysis also let us discuss the 
charge transfer (CT) from the T to S atoms, which also increases with the applied strain. 
Therefore, the T-S bonds become more ionic in comparison with the equilibrium structures. 
In addition, the differences in the atomic charges of both metals suggest that M0S2 is more 
polarized compared with WS2. 

The geometry changes due to the applied deformations are elastic and the materials 
should return to their equilibrium structures once the strain is released. Compliance with 
Hooke's law is confirmed by our calculated stress-strain relation curves, which show linear 
scaling within the range of applied e. The stress-strain curves for the uni-axial deformations 
are shown in Fig. 3. The observed linear relation means that we are working in the truly 
elastic region, and that the TMDs are very ductile materials. Fitting the linear function to 
the stress-strain curve gives the corresponding elastic modulus. For the M0S2 monolayer, 
this elastic modulus is 208 GPa and 182 GPa for the armchair and zigzag directions, respec- 
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e = AL/L Q e = AL7L Q e = AL/L Q 

Figure 2: (Online color) Calculated T-S bond lengths (a), S-T-S bond angles (b), and the 
T atomic charges (c) of M0S2 (top) and WS2 (bottom) monolayers under tensile strain and 
compression. Full symbols denote the variables along the deformation axis, empty symbols 
correspond to the variables along the fixed axis. 



6 




e = Al_/I_ e = AL/L 

Figure 3: (Online color) Calculated strain-stress relation of M0S2 (a) and WS2 (b) mono- 
layers under uni-axial tensile strain. The Young's moduli for the deformations along the a 
and b lattice vectors are indicated by Y £a and Y ez for the armchair and zigzag edge types, 
respectively. 



tively. Our results are in very good agreement with the recent work of Li et al., 31 who have 
reported 200.3±3.7 GPa and 197.9±4.3 GPa, respectively using DFT/PW91 level of theory 
as implemented in PWscf code. We note that our elastic modulus is 22 GPa smaller than 
the Young's modulus reported by Lorenz et al. 32 for the M0S2 nanotubes calculated at the 
DFTB level. This deviation can be related to the differences in the optimized monolayer 
thickness and the quantum-mechanical method used for the computations. The stress-strain 
relation and the elastic modulus obtained for WS2 is nearly the same for the zigzag and 
armchair directions (217 vs. 211 GPa). 

Fig. 4 shows the band gap evolution with the applied 2D-isotropic strain for M0S2 and 
WS2 hexagonal monolayers. In equilibrium, the single TMD layers are direct band gap semi- 
conductors with the transition at the K point. Under tensile strain, the band gap decreases 
linearly and the system exhibits an indirect transition from the V to the K point. Eventu- 
ally, the band gap vanishes for elongations of about 11%, and consequently the monolayers 
become metallic with bands crossing the Fermi level at K (conduction band) and Y (va- 
lence band) points. At this elongation the T-S bonds are stretched, but not yet broken (cf. 
Fig. 2). Similar behavior is observed when a monolayer is exposed to compressive strain, 
though the influence is smaller than that of the tensile strain. Uni-axial tensile strain results 
in a similar band gap lowering with eventual transition to the metallic state. In this case, 
however, for a given amount of elongation or compression, the changes are smaller than those 
observed for 2D-isotropic deformations. These results are in a very good agreement with 
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Figure 4: Band gap evolution of M0S2 (a) and WS2 (b) monolayers under isotropic tensile 
strain and compression. The fc-points, at which the top of valence band and the bottom 
of conduction band occur, are indicated. The equilibrium structures are marked with red 
circles. 



the works of Johari et al., 13 Yue et al., 14 and Peelaers et al. 33 The influence of mechanical 
deformations on the electronic structure is mostly visible at the K point, where both the 
bottom of conduction and the top of valence bands are significantly lowered compared with 
the equilibrium situation. At the same time, the top of valence band at the Y point moves 
upwards (see Fig. 5). This behavior can be understood as the orbital overlap between the 
Mo and S atoms decreases due to the change in the bond distances. 14 

The changes in the electronic structures under mechanical deformations can enhance or 
quench the electron conductivity through semiconducting TMD monolayers. Fig. 5 shows the 
comparison of electronic structures and the corresponding electron conductivities of M0S2 
and WS2 monolayers. Three examples are plotted, namely direct band gap in the equilibrium 
structure, indirect band gap (4.7% elongation), and the metallicity (12.6% elongation). The 
2D-isotropic strain (e) gives rise to reduction of the band gap and significant increase in the 
electron conductivity at the Fermi level. 

Since at approximately 11% deformation band gaps of M0S2 and WS2 vanish, the trans- 
port channels are completely open and electrical conductivity is enabled. This can be ex- 
plained by the deformation of the Brillouin zone, where the applied strain uniformly changes 
the primitive cell, bands cross the Fermi level, and consequently the character of a mono- 
layer changes from the semiconducting to metallic. The electronic transport in the case of 
isotropic strain/compression does not depend on the direction of the current and it is almost 
unchanged for the zigzag (Q z ) and the armchair (Q a ) directions. 

Fig. 6 shows the electron conductivity evolution of M0S2 monolayers with uni-axial strain 
applied in the parallel (e||) and perpendicular (e±) directions to the transport axis. For 
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Figure 5: (Online color) Band structure, the corresponding density of states, and the electron 
conductivity of M0S2 (a) WS2 (b), and NbS2 (c) monolayers under isotropic tensile strain. 
Q a and Q z denote electron conductivity along the armchair and zigzag direction, respectively. 



elongations as large as 13%, the monolayers still show semiconducting behavior. This case 
corresponds to the stretching of TMD nanotubes with very large diameters and is consistent 
with our results on the mechanical deformation of these materials 34 and the recent results 
of Lu et al. 17 Nevertheless, e± shows a very useful way of tuning the band gaps of TMD 
monolayers in a wide range without changing the conductivity significantly (see Fig. 6 a). 
In contrast, eii behaves differently and a decrease in the band gap result also in the reduced 
conductances for the electron energies lower than the Fermi energy. At the energies higher 
than the Fermi level, £u enhances the conductance dramatically (see Fig. 6 b). These results 
highlight the importance of the strain direction for nanoelectromechanical applications. 

For comparison, we have also calculated the electronic transport of metallic NbS2 mono- 
layers under tensile strain. At the Fermi level, the electron conductivity is present indepen- 
dent of the mechanical deformation, however, the intensity decreases significantly for both 
isotropic and uni-axial tensile strain (cf. Figs. 5 and 6). 

4 Conclusion 

Using non-equilibrium Green's functions (NEGF) 20 method combined with the Landauer- 
Biittiker approach we have calculated quantum transport of the mechanically deformed M0S2 
and WS2 monolayers. Under 2D-isotropic strain, the TMD monolayers change the semicon- 
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Figure 6: (Online color) Electron conductivity of M0S2 (a) and NbS2 (b) monolayers under 
uni-axial tensile strain. 



ducting direct gap character and become metallic for the deformations as large as 11%, what 
results in electrical conductivity at the Fermi level. Uni-axial tensile strain of monolayers 
has less influence on the electronic structure for a given elongation. This corresponds also 
to the deformations of large-diameter TMD nanotubes with zigzag and armchair chiralities 
for the e z and e a , respectively. The semiconductor-conductor transition is reflected in the 
electron transport, our calculations show similar values for metallic NbS2 as for stretched 
M0S2 and WS2. We expect that our calculations for monolayers are directly transferable 
to TMD nanotubes, in particular for the experimentally available ones with rather large 
diameters. A study on TMD nanotubes is in progress and will be a subject of a forthcoming 
communication. 34 

Generally, mechanical deformation is a possible way of tuning the band gaps and trans- 
port properties of TMD monolayers due to the metal-semiconductor transition. Therefore, 
strain effects in TMDs can opens new prospects for their applications in optoelectronics or 
nanosensor device engineering. 
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Support Infromation 

Band structures and partial densities of states of TS2 (T = Mo, W) monolayers under 
mechanical deformations calculated at the DFT/PBE level (see Figs. 1 and 2). 
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Figure 1: (Online color) (a) Band structures of equilibrium and deformed M0S2 monolayers 
calculated at the DFT/PBE level. The horizontal dashed lines indicate the Fermi level. 
The arrows indicate the fundamental band gap (direct or indirect) for a given system. The 
top of valence band (blue) and the bottom of conduction band (green) are highlighted, (b) 
Corresponding partial densities of states of M0S2 monolayer with projections on the Mo and 
S atoms: the contributions from Ad and 3p orbitals of Mo and S, respectively. 
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Figure 2: (Online color) (a) Band structures of equilibrium and deformed WS2 monolayers 
calculated at the DFT/PBE level. The horizontal dashed lines indicate the Fermi level. 
The arrows indicate the fundamental band gap (direct or indirect) for a given system. The 
top of valence band (blue) and the bottom of conduction band (green) are highlighted, (b) 
Corresponding partial densities of states of WS2 monolayer with projections on the Mo and 
S atoms: the contributions from 5d and 3p orbitals of W and S, respectively. 
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